Polarization and decoherence in a two-component Bose-Einstein Condensate 
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We theoretically investigate polarization properties of a two-component Bose-Einstein condensate 
(BEC) and influence of decoherence induced by environment on BEC polarization through intro- 
ducing four BEC Stokes operators which are quantum analog of the classical Stokes parameters for 
a light field. BEC polarization states can be geometrically described by a Poincare sphere defined 
by expectation values of BEC Stokes operators. Without decoherence, it is shown that nonlin- 
ear inter-atomic interactions in the BEC induce periodic polarization oscillations whose periods 
depend on the difference between self-interaction in each component and inter-component interac- 
tion strengths. In particular, when inter-atomic nonlinear self-interaction in each BEC component 
equals inter-component nonlinear interaction. Stokes vector associated with Stokes operators pre- 
cesses around a fixed axis in the dynamic evolution of the BEC. The value of the processing frequency 
is determined by the strength of the linear coupling between two components of the BEC. When 
decoherence is involved, we find each component of the Stokes vector decays which implies that 
decoherence depolarizes the BEC. 
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I. INTRODUCTION 

The experimental realization of Bose-Einstein conden- 
sation ITj and an atom laser Oj have sparked many theo- 
retical and experimental studies of coherent atomic mat- 
ter. It is well understood that matter- wave analog of 
nonlinear optical interactions is provided by inter-atomic 
collisions. In particular, in the s-wave scattering regime, 
two-body collisions are mathematically equivalent to a 
optical Kerr interaction. Hence, many of the concepts 
of first developed in optics can readily be extended to 
Bose-Einstein condensates (BEC). Recently, much atten- 
tion has been paid to these lines, including the study of 
matter- wave solitons 0, phase conjugation Q, four- wave 
mixing (sj, an atom laser atom holography 0, etc. 
Especially, pulsed and continuous-wave atom lasers 
and the four-wave mixing with matter waves have 
been demonstrated experimentally. The experimental 
verification of the other predictions is also underway. 

As is well known that polarization is an important 
property that is common to all types of vector waves. 
Light waves possess this property and so do elastic and 
spin waves in solids, for example. In this paper we study 
BEC polarization that arises for matter waves by con- 
sidering the two components of a two-component BEC 
as two polarizations. The motivation for present study 
is twofold. First, the BEC polarization is interesting in 
itself and deserves further examination. From a basic 
point of view, it significantly extends the already well- 
established range of analogies between light and mat- 
ter waves. Second, more importantly perhaps, it may 
also have useful practical applications. We shall inves- 
tigate nonlinear polarization dynamics of the BEC and 
study influence of decoherence induced by noncondensate 
atoms. 



To begin we consider a two-component BEC described 
by the second quantized Hamiltonian 



Hi+ H2+ Hni + Hi , 



H 

= y"rfxv31(x)[-|^v2 + y,(x) 

+C/,^|(x)V-,(x)]V3.(x),(^ = 1,2), 

Hnl = Ui2 J dxV3t(x)^^(x)^i(x)^2(x), 

Hi = A f dx[7Al(x)V'2(x)+^i(x)7At(x)]. 



(1) 

(2) 
(3) 

(4) 



Here i = 1,2, ^pi{x) and ?/'J(x) are the atomic field op- 
erators which annihilate and create atoms at position 
X, respectively. Hi and H2 describe the evolution of 
each component in the absence of intercommunion in- 
teractions. Hnl is a nonlinear interaction term, it de- 
scribes inter-component collisions. Hi is a linear inter- 
action term, it describes an atomic Josephson tunnelling 
process. I^(x)(i = 1,2) are trapping potentials. Inter- 
actions between atoms are described by a nonlinear self- 
interaction term Ui — ATrh?af'^ /m and a term that corre- 
sponds the nonlinear interaction between different com- 
ponent U12 = A-nh^a']^/m, where af^ is s-wave scattering 
lengths of component i and 2 that between component 
1 and 2. For simplicity, throughout this paper we set 
h=l, and assume that = al"^ = a^'^, Vi(x) = V2(x). 

For weak many-body interactions, i.e., for small num- 
ber of condensate atoms, it has been well known that the 
Hamiltonian (1) can reduce to a two- mode Hamiltonian 



H 



ujf){a\ai + 02^2) + q{a\ + a\ 0,2) 
+g{a\a2 + a\ai) + 2xa\aia\a2, 



(5) 
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where {af) is correspondent-mode annihilation (cre- 
ation) operators with [a^, a|] = 1, and q, x and g are cou- 
pling constants which characterize the strength of non- 
hnear and hnear interactions, respectively. 

In general, the Hamiltonian ©can not be exactly 
solved except the two limiting cases {i) g — 0, q ^ 
and/or X 0; (ii) 5 and g = x = 0. However, 
for the general case of g ^ 0, q ^ and x 7^ 0; an ap- 
proximate analytic solution of the Hamiltonian jSJ can 
be obtained by introducing a pair of bosonic operators: 



-igt. 



ai 



+ a2)/V2 and A2 



0.2 



)/V2, 



which satisfy the usual bosonic commutation relation: 
= 5ij. Assume that the BEC is initially in the 
coherent state defined by 



|ai,a2) = exp [-/racl2(|aip + jaip)] ^ 



tj_u.2 



n,m— U ^ 



\n, m) 



(6) 



Where ai and a2 are two arbitrary complex numbers. 
Under the rotating wave approximation, the energy and 
wave function of the system in the (Ai, A2) representa- 
tion with the basis {\n, m)} can be given by 



Point representing the 
polarization state 




LCP I 



E{n,m) — Lu{n + m) + g(n — m) + -{3q + 2x) 

X (n + m)^ — ^ (n — m)^ — x^"^i C'') 



FIG. 1: Representation of the BEC polarization state on the 
Poincare sphere. 
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(8) 



where n and to take nonnegative integers, a; = — (x ~ 
g)/2, = laip-f |Q;2p = |uip-h|M2p is the total number 
of atoms with two new parameters defined by 

Ui — —{ai+a2), U2 — —^{ai — a2). (9) 



V2' 



V2' 



II. BEC POLARIZATION STATES 



can be completely described by the four Stokes param- 
eters which can be written as a string: {So, Si, S2, S3}. 
The first parameter expresses the total number of atoms 
in the BEC. The remaining three parameters describes 
the polarization state. 

For a completely unpolarized BEC only the parameter 
is nonzero: 



5^ = 5o{l;0;0;0}. 



(11) 



The BEC is called completely polarized if the Stokes 
parameters satisfy the following relation: 



5'o — {Si + S2 + s'^y 



(12) 



In order to discuss the polarization properties of the 
two-component BEC, it is convenient to define the Her- 
mitian Stokes operators ^3 as 

•S'o = (0^01+0^02), = (a+oi - 0^02), 

S2 = (a^ a2 + ai) , 53 = 1(0^01-0+02), (10) 

which satisfy the commutation relations: [Sj,Sk] — 
2iejkiSi and [5*^,50] = 0. The noncommutability of the 
Stokes operators Si, S2, and 6*3 precludes the simultane- 
ous measurement of the physical quantities represented 
by them. 

These Stokes operators are quantum analog of the clas- 
sical Stokes parameters for a light field, i.e., 5^ = (5*^), 
where fi = 0, ...,3. Any polarization state of the BEC 



A partially polarized BEC is a mixture of a completely 
unpolarized BEC and a completely polarized BEC. All 
of the four Stokes parameters are independent, however, 
they satisfy the following inequality 5*0 > {Sf -I- 5*1 -I- 
S^y^"^. And the degree of polarization can be measured 

by 



P 



VS't + Si + S: 



So 



(13) 



For a partially polarized BEC < P < 1, a completely 
polarized BEC has P = 1, while a completely unpolarized 
BEC has P = 0. 

BEC polarization states can be described by a Poincare 
sphere plotted in Fig. 1. A particular polarization 
state of the BEC may be represented as a point in a 
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three-dimensional Stokes space by introducing the Stokes 
vector S with Cartesian coordinates Sj, S2 and S3. 
The point hes on a Poincare sphere jlflj with radius 
{Sf + S2 + Sfy/^. If the Stokes parameters are known, 
the polarization azimuth and the ellipticity angle rj 
defined by the following expressions, respectively. 



tan 29 
sin 2ri 



S2 
Si' 
S3 
So' 



(14) 
(15) 



If the total number of the atoms, and therefore Sq, the 
degree of polarization P, the polarization azimuth O and 
the ellipticity 77 of the polarization ellipse are known, one 
may return the Stokes parameters 5i_3: 



51 = PS'oCOs(2r/)cos(2e), 

52 = PS'oCOs(2r/)sin(2e), 
^3 = P^osin(2?7). 



(16) 



If the BEC is totally polarized, the radius of the 
Poincare sphere equals So- The 'north' and 'south' poles 
of the sphere correspond to right and left circular polar- 
izations, respectively. Points on the 'equator' represent 
linear polarization with 'longitude' being twice the po- 
larization azimuth @. The 'latitude' of a point will give 
twice the angle of ellipticity r]. Antipodes, i.e., points 
at opposite ends of a diameter of the sphere, are known 
as orthogonal polarizations. If the BEC is partially po- 
larized the ratio between the Poincare sphere radius and 
the Stokes parameter Sq gives the degree of polarization 

m- 

We now examine the dynamical evolution of polariza- 
tion of the BEC. From Eqs.Q and © we can calculate 
the Stokes parameters. The results are 

^0 = N, (17) 
= 2\ui\\u2\cos[4:gt + e{t)] 

xexp[~-2Nsm^{'^-^)tl (18) 

5*2 = 2|ai||a2| cos(i^ai - l^aj, (19) 
^3 = -2\ui\\u2\sm[4^gt + e{t)] 

xexp[~2Nsm^{'^^)t]. (20) 



where we have set Ui = \ui\e^''^"i and the symbol: 

0{t) = {^,,,~^,,,) + {\u2\^-\ui\^)sin{q-x)t, (21) 

Assume that the system is initially in a state jai, > 
of ai being real and 012 = 0, which implies that the BEC 
is initially in a linear horizontal polarization state with 
the Stokes parameters given by 5*^ = S'ojl; 1; 0; 0}. From 
Eas. H17|l to H19() . one can obtain some interesting polar- 
ization states of the BEC at some specific times. For 
instance, at times = (2fc-t- l)7r/ {8g) with k being an in- 
teger, the BEC is in a left-hand circular polarization state 




FIG. 2: Precession of the Stokes vector on the Poincare sphere 
in the BEC evolution when q — x 8-nd without decoherence. 



with the Stokes parameters being S"^ = S'o{l; 0; 0; — 1}. 
When tk = {2k — l)7r/(8.g) with k being an integer, the 
BEC is in a right-hand circular polarization state with 
the Stokes parameters being 5^ = So{l; 0; 0; 1}. At times 
tk — kn / (2g) with k being an integer, the BEC returns 
the initial polarization state, a linear horizontal polariza- 
tion state. When tk = {2k + l)TT/{Ag) with k being an 
integer, one may find that the BEC is in a linearly ver- 
tical polarization state with the Stokes parameters given 
by5^ = 5'o{l;-l;0;0}. 

From Eqs. (11) to (14), we can obtain the degree of 
polarization 



1 



Pit) = j^{{\u2\'^\ui\^r+Mui\'\u2\' 



X exp 



1/2 



(22) 



From above equation we see that: (i) the BEC system 
exhibits periodic polarization oscillations with a period 
T = An/lq — x\ which depends on the difference between 
the self-interaction and inter-component interaction, (ii) 
the degree of polarization depends only nonlinear inter- 
actions and the initial state of the BEC but not linear in- 
teraction due to the independence of P{t) upon the linear 
couphng. (iii) When q = x, we obtain P{t) = P(0) = 1, 
which means that the degree of polarization remains un- 
changed. It is interesting to note that when q = Xy from 
Eqs. H17|) to (|19|l . we can find that the evolution equation 
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of the Stokes vector takes the following form 



we can get the Stokes parameters 



dS 



— = -rj X S, 



(23) 



where fl — 4:ge2 with 62 being the unit vector of the 
S'2 axis. Geometrically, Ea. H23|l indicates that the Stokes 
vector S precesses around the S'2 axis in the dynamic evo- 
lution of the BEG. The value of the processing frequency 
il is determined by the strength of the BEG linear cou- 
pling. In Fig. 2 we plot the precession of the Stokes vector 
in the BEG evolution. 



III. INFLUENCE OF DECOHERENCE ON 
POLARIZATION 



We now consider the effect of the decoherence on 
the BEG polarization. As is well known, any system 
inevitably interacts with its environment. In fact, in 
current experiments on BEGs of dilute alkali atomic 
gases condensate atoms continuously interact with non- 
condensate atoms (reservoir). Interaction between a 
BEG and environment causes decoherence [Hill. We 
use a reservoir consisting of an infinite set of harmonic 
oscillators to model environment of condensate atoms in 
a trap, and assume the total Hamiltonian to be 



H 



F{{d})Y,c,ibi + k) 



(24) 



where the second term is the Hamiltonian of the reser- 
voir. The last term in Ea. (|24|) is a renormalization term 
[1^ . The third term in Eq. H24() represents the interaction 
between the system and the reservoir with a coupling 
constant Cfc, where {O} is a set of linear operators of the 
system or their linear combinations in the same picture 
as that of H, F{{0}) is an operator function of {O}. In 
order to enable what the interaction between the system 
and environment describes is decoherence not dissipation 
[15|. we require that the linear operator O commutes 
with the Hamiltonian of the system H. The concrete 
form of the function F{{0}), which may be considered 
as an experimentally determined quantity. Therefore, the 
decohering interaction in l|24|) can not only describe de- 
coherence caused by the effect of elastic collisions be- 
tween condensate and non-condensate atoms, but also 
simulate decoherence caused by other decohering sources 
[T^ through properly choosing the operator function of 
the system F({6}). 

From Hamiltonian (|24l) through lengthy calculations 



So{t) = N, (25) 
Si{t) = 2^sin'(/)J„„(t)exp[-7(™_„)(™+i^„_i)(t)];,26) 

S2it) = ^cos"(/)^„(i)exp[-7(„,„)(„,„)(i)], (27) 

m,n 

Ssit) = ^cos>;„(t)exp[-7(„^„)(^+i_„_i)(i)] 

m,n 

- sin" (f>'^^{t) exp[-7(„^„)(„^„)(i)]. (28) 

rn.n 

where we have introduced symbols 



Sin''/''„„(t) = V«^™+l)IP(m,n)(m+l,n-l)(0)l 

X sin[2.gi + <l>(m,n)(m+im-i){t)], (29) 
cos' (/)mn(0 = y/n{m + l)|p(,„,„)(,„+i,„-i)(0)| 

X C0s[2gt + (t){m,n){m+l.n-l)it)], (30) 
sin'^mnW = (W- »^)|P(m,n)(m,n)(0)| 

X Sm[2gt + (l)(m,n)(m,n){t)], (31) 
cos" (j)'^n{t) = ("I - ?^)|P(m,n)(m,n)(0)| 

X COs[2gt + (j)(^,n,n){m.n){t)], (32) 

where the damping factor and the phase shift induced by 
interaction between the BEG and reservoir are given by 

l{m',n')(m,n){t) = v'i{m' , u' \ m, n)Q2{t) , (33) 
<l>{m',n')(m,n){t) = v+{m' , Ti' ; TJi, n)v-{m' , Ti' ; m, n)Q i{t) 

+&im\n'){m,n), (34) 

Here we have introduced the notations 
v±{n,m;n',m') = F({0(n, m)}) ± F({0(n', m')}035) 

(0) = \P{,n',n' )(m,n) (0) 

x|exp[-i6l(,„/_„/)(,„^„)]. (36) 



where the two reservoir-dependent functions are given by 

Qi{t) = 



dwJ(w)^--^sin(wt), (37) 
^ 



Q2{t) = 2 / 
Jo 



Lot 



X sin^ — coth — , 



(38) 



where J{lo) is the spectral density of the reservoir, c(w) is 
the continuum expression for Cfc, and (3 = l/ksT with ks 
and T being the Boltzmann constant and temperature, 
respectively. 

It is well known that decoherence corresponds to the 
decay of off-diagonal elements of the reduced density ma- 
trix of a quantum system. For the case under our con- 
sideration, the degree of decoherence is determined by 
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the decaying factor given by Ea. (|33|l . It is interest- 
ing to note that if we choose a proper operator func- 
tion F{{d}) to make F{{0{m' ,n')}) = F{{0{m,n)}) 
for (rn',n') ^ {m,n), then we find that p(m' ,n'){m,n){i) — 
P(m' ,n')(m,n){^) which indicates that the decoherence-free 
evolution of the BEC is reaHzed, and the BEC preserves 
its initial polarization state. From Eq. H26(l into H28(l we 
can immediately draw one important qualitative conclu- 
sion: since Q2{t) is positive, the existence of the decoher- 
ence is always to tend to make Stokes parameters decay. 
Therefore, decoherence always depolarizes the BEC. 

From Eqs.lUnjl to , ^ and ^ we see that aU 
necessary information about the effects of the reservoir 
on the polarization is contained in the spectral density 
of the reservoir. To proceed further let us now specialize 
to the Ohmic case |[17j with the spectral distribution 

where uJc is the high frequency cut-off, k is a posi- 
tive characteristic parameter of the reservoir. With this 
choice, when uJct S> 1, at zero temperature we have 
Q2{t) = K\n{ujct), so that the components of the Stokes 
vector are given by 

S,it) = 2^sin>:„„(i)M)-'''^-("'"^'"+^^"-^\(40) 

7n.n 

S^it) = ^cos"C„W(^ct)-«"-(™'"^"'"\ (41) 
+ cos' 0:„„(t)(wct)~""-^"'"'™+''""'^(42) 

m.n 

which indicate that Stokes parameters decay according 
to the 'power law'. Therefore, geometrically decoherence 
makes the Stokes vector contract towards the center of 
the Poincare sphere according to the 'power law' with the 
time evolution. 

At finite temperature, we have Q2(0 ^ + 
so that 

S^{t) = 2^sin'0:„„(t)(^)---("."^"+i^"-i) 

m,n 

X exp[— K7r/3~"'^w^(m, n;m + l,n — l)t], (43) 
S2{t) = ^cos"C„(t)(^)"--(™'"^™.«) 

771 . n 

X e:>cp[—K.7Tp~^v^{m^ n; m, n)t], (44) 
S,it) = ^cos>'_(t)(§:^)— 

^ — ' ZTT 

m,n 

X exp[—KTr(3~^v'^{m, n; m + 1, 71 — l)t], 
-^sin"C„(i)(^)^'''^-<'"'"^"'"^ 

^ — ' ZTT 

m,n 

X exp[— K7r/3~"'"u^(m, n; to, n)t], (45) 



which indicate that at finite temperature Stokes param- 
eter decay according to the 'exponential law'. Geometri- 
cally this means that decoherence makes the Stokes vec- 
tor contracts toward the center of the Poincare sphere 
according to the 'exponential law' with the time evolu- 
tion. Therefore, decoherence depolarizes the BEC. 

From Eqs. (gOJ to we see that a BEC initially 
in a completely polarized state (P = 1) evolves into a 
completely unpolarized state (P = 0) when the time ap- 
proaches the infinity due to the influence of decoherence. 



IV. CONCLUDING REMARKS 



In conclusion, we have examined the nonlinear polar- 
ization dynamics of a two-component BEC and influence 
of decoherence induced by environment on BEC polar- 
ization through introducing four BEC Stokes operators. 
We have shown that different polarization states of the 
BEC may be geometrically represented on a Poincare 
sphere. Without decoherence, it has been found that 
nonlinear inter-atomic interactions in the BEC induce 
periodic polarization oscillations whose periods depend 
on the nonlinear interaction strengths. When nonlinear 
self-interaction (g) equals inter-component nonlinear in- 
teraction (x) , the BEC initially in a completely polarized 
state remains a completely polarized state, however, the 
Stokes vector precesses around the S2 axis with a precess- 
ing frequency determined by the linear coupling strength. 
When decoherence is taken into account, we find each 
component of the Stokes vector decays. This indicates 
that decoherence depolarizes the BEC. In particular. We 
have obtained analytical expressions of the Stokes pa- 
rameters and found that for the reservoir-spectral den- 
sity of the Ohmic case, the stokes parameters decay by 
the 'power law' at zero temperature, and the 'exponential 
law' at finite temperature, respectively. 
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